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The detection and the measurement of the production cross section of a light Higgs boson at the 
CERN Large Hadron Collider demand the accurate prediction of the background distributions of 
photon pairs. To improve this theoretical prediction, we present the soft-gluon resummed calcu- 
lation of the pp — > 77 A cross section, including the exact one loop gg — > 7717 contribution. By 
incorporating the known fixed order results and the leading terms in the higher order corrections, 
the resummed cross section provides a reliable prediction for the inclusive diphoton invariant mass 
and transverse momentum distributions. Given our results, we propose the search for the Higgs 
boson in the inclusive diphoton mode with a cut on the transverse momentum of the photon pair, 
without the requirement of an additional tagged jet. 



PACS number(s): 12.38.Cy, 13.85. Qk. 



I. INTRODUCTION 

The direct search for the Standard Model (SM) Higgs 
boson at the CERN LEP collider constrains its mass ijih 
to be greater than y/S — m z ~ 108 GeV Q. Meanwhile, 
recent electroweak global fits , and the measurements 
of the W ± boson and top quark masses ||] suggest that 
the SM Higgs mass is lower than about 260 GeV. If run is 
less than twice the Z° boson mass, as the electroweak fits 
hint, then at the CERN Large Hadron Collider (LHC) 
the most promising detection modes of the SM Higgs 
boson will be pp — > HX — > 77 AT Q and the associated 
production pp — > H jet X — > 77jetAT 0]. According to 
earlier studies, a statistical significance on the order of 5- 
10 can be reached for both of these signals, actual values 
depending on luminosity and background estimates. In 
Ref. H it was also found that in order to optimize the 
significance it is necessary to impose a 30 GeV cut on the 
transverse momentum of the jet, or equivalently (at NLO 
precision), on the Qt of the photon pair. With this cut in 
place extraction of the signal in the Higgs plus jet mode 
requires the precise knowledge of both the signal and 
background distributions in the mid- to high-Qy region. 

The precise determination of signal and background 
rates of the inclusive diphoton process demands the cal- 
culation of the large QCD corrections. Higher order cor- 
rections, both to the signal H and the background |t]], in- 
crease the rate significantly, by a factor of about 2. In the 
case of the background this large increase is mostly due 
to the fact that the diphoton cross section receives a large 
contribution from the formally higher order gg — ► 77 par- 
tonic subprocess, which is of the same order of magnitude 
as the qq + qg — > 77 AT subprocesses [|| . Since the lowest 
order gg — > 77 subprocess proceeds through a box dia- 
gram, the calculation of yet higher order corrections to 
this process is complicated. 
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A reliable calculation of the transverse momentum dis- 
tribution of the Higgs boson or the background photon 
pair also requires the resummation of the potentially 
large logarithmic contributions of the type h\{Q/Qr) 
(where Q is the invariant mass of the pair), arising as 
a result of multiple soft-gluon emission. Using the soft- 
gluon resummation technique, the low- to mid-Qr region 
can be predicted, and the resummed calculation can be 
matched to the fixed order, giving a reliable prediction 
in the whole Qt range (|||. As an added bonus the re- 
summed cross section also exhibits reduced scale depen- 
dence since it includes the most important higher order 
contributions. It also gives a hint of the size of the yet 
higher order corrections. 

The effects of the multiple soft-gluon radiation on the 
transverse momentum distribution of Higgs bosons were 
discussed in a recent paper In the present work, we 
analyze the same effects on the transverse momentum of 
the background photonpair, extending results previously 
published in Refs. [§,[l(J. When Qt is integrated, the 
resummed calculation of Ref. || gives the rates for the 
QQ + Q9 ~ * 77 A" subprocesses at the 0(a 2 as) precision. 
That calculation also includes the photon fragmentation 
contribution, and approximates the 0(a 2 a^) gg + qg — > 
77AT rate. In this work, we include the exact one loop 
0(a 2 a^) gg — > 7717 calculation to improve the diphoton 
background prediction in the high Qt region. The above 
fixed order and resummed contributions are implemented 
in the ResBos IQJ Monte Carlo event generator which 
was used to produce our numerical results. 



II. ANALYTICAL RESULTS 

The one loop gg — > 7717 amplitude can easily be de- 
rived from the one loop five-gluon (5<?) amplitudes of 
Ref. |Lq] for the case when the particles in the loop are 



fermions in the fundamental representation, by replacing 
two of the gluon vertices with photon vertices. Since the 
5g amplitude is explicitly given in a color-decomposed 
form, it is possible to replace the SU(3) generators and 
strong couplings (gs) of two of the gluons with U(l) gen- 
erators and the electromagnetic couplings (e) of photons. 
The final expression for the square of the three-gluon- 
two-photon (3^27) amplitude is 



\A(g l9 2 S37475)r = 8(eQ l ) 4 5s^c(^ - 1) x 



E 



helicities 



E ^5*1 ( CT l)°'2,O'3,0'4,0'5) 



cr£GOP\ 



(1) 



where cr, is shorthand for the 4-momenta and helicities, 
{Pm, ^cn}, of the gluons 1,2,3 and the photons 4,5. The 
charge of the quarks in the loop is given by Qi in the units 
of the charge of the positron, and Nc = 3 is the number 
of colors in QCD. In Eq.(Q) COP^ 123 ' denotes the subset 
of permutations of (1,2,3,4) that leave the ordering of 
(1,2,3) unchanged up to a cyclic permutation 
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A 5:l ( CT 1 ^2, T3, 0"4, 0- 5 ) 



Ag{ 9] (l, 2, 3,4,5) + Ag{ a] (1,2, 4,3,5) + 
A${ a] (1,4, 2, 3, 5) + 4-f ] (4,l, 2, 3, 5) + 
4^3, 1,2, 4, 5) + 4^ 2] (3, 1,4, 2, 5) + 
4/ 2l (3,4,l,2,5) + 4( 2l (4,3,l,2,5) + 
4^ 21 (2,3,1,4,5) + 4^ 2] (2,3,4,1,5) + 
Agf 1 (2,4,3,l,5) + 4^ 2] (4,2,3,1,5). 

[1/2], 



(2) 



The partial amplitudes A^.[ 1 (1, 2, 3, 4, 5) for the various 
helicities of the external gluons and photons are given in 
Ref. [@. 

The extension of the Collins-Soper-Sterman soft-gluon 
resummation formalism Jl3| to diphoton production was 
discussed in Ref. [[l0| . We follow that work when calculat- 
ing the 0(as) fixed-order corrections to the qq — > 77 sub- 
process, including the fragmentation contributions, and 
resumming the contributions of the qq — > 773, qg — > 779 
and qg — > 779 subprocesses. In addition, Ref. |l(J used 
an approximate expression for the cross-section of the 
99 ~ * 775 subprocess, that accounted for the emission 
of soft and/or collinear gluons off the initial gluons, and 
was valid in the limit of small non-zero Qt- 

Using the exact formula (Q) for the matrix element of 
the gg — > 7717 subprocess, we can now also analyze the 
high-Qy region. We match the above fixed-order result 
with the resummed results of Ref. Jl0| , following the pro- 
cedure of Ref. jy]], to obtain a prediction in the whole 
range of Qt- In the resummation of the gg + qg — > 773 
process, we use the A^ and A^ coefficients of Ref. [n4j, 
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FIG. 1. The resummed part of invariant mass and trans- 
verse momentum distributions of photon pairs at the LHC, 
calculated for the qq + qg — > J"/X subprocess using ResBos. 
The solid curves are calculated using the exact Wilson coeffi- 
cient C (1) . For the dashed curves, C' 1 ' is set identical to that 
of the Drell-Yan process. For the dotted curves, is set to 
zero. 



since these coefficients depend only on the initial par- 
tons and are independent of the partonic process itself. 
The coefficient for the gg — ► 7717 process is the 

same as the one for gg — > Hg 7 also given in Ref. p4[ . 
The only unknown part of the resummed cross-section 
for this subprocess is the coefficient C^> , which would 
require the computation of two-loop virtual diagrams. 
Following Ref. |lC| ], we approximate in the process 
99 + 99 ~* 775 by the expression for in the process 
gg — > H g in the limit of a large top quark mass. This ap- 
proximation is reasonable when both photons have large 
transverse momenta, because then the real and virtual 
corrections are dominated by soft and collinear radiation 
off external lines, which is identical to that in the Higgs 
production process. 

To support this statement, we examine the behav- 
ior of the analogous approximation for the case of the 
qq, qg, qg — > "f"fX process. We calculate the resummed 
part of this cross section in three different ways: with 
an exact C«, with an approximate C' 1 ', and without 
the C' 1 ' coefficient. The approximate coefficient 
is identical to that of the Drell-Yan process, and omits 
terms coming from the additional virtual corrections to 
the hard process (c.f. Eqs. (8) and (11) of @). The 
results are displayed in Fig. [y. The figure shows that 
the approximation (dashed curve) is better toward the 
high invariant mass region, which is expected, since Q 
is correlated with the transverse momentum of the indi- 
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vidual photon in the central rapidity region. According 
to the invariance mass plot, the approximation slightly 
overestimates the rate for low Q's (by a few percent), and 
the (dotted) curve without deviates from the exact 
(solid) curve by about 20 percent, which is the size of the 
NLO corrections to this process. In conclusion, the ap- 
proximate CW coefficient (borrowed from the Drell-Yan 
process) captures most of the NLO corrections for the 
qq, qg, qg — > 77 A process, and is a good ansatz for both 
Q and Qt, when compared with the exact calculation. 
In the absence of the exact calculation for the co- 
efficient, which requires the knowledge of the exact two 
loop virtual corrections to the gg + qg — > 77 process, 
we propose to use the approximate C^ 1 ' coefficient, bor- 
rowed from the gg — > HX process, to estimate the effect 
of the higher order corrections to the distributions of the 
photon pairs. We expect that this should give a bet- 
ter estimate of the event rates than using only the 
coefficient in our calculation. 

Finally, the small qg component is still approximated 
as in Ref . |l(| , noting that it is highly suppressed due to 
the large difference between the gluon and quark lumi- 
nosities in the probed region of momentum fraction. 



III. NUMERICAL RESULTS 

The above described analytic results are implemented 
in the ResBos Monte Carlo event generator [^T). In the 
numerical calculations we use the center-of-mass energy 
14 TeV and the following electroweak parameters |H| : 



G F = 1.16639 x 10" 
m w = 80.41 GeV, 



GeV , m z = 91.187 GeV, 

a(mz)= . 

v ; 128.88 



We use the NLO expressions for the running electromag- 
netic and strong couplings a{p) and as(p), as well as the 
NLO parton distribution function set CTEQ4M. We set 
the renormalization scale equal to the factorization scale: 
fj,R = fiF = Q. For the choice of the non-perturbative 
function of the resummation for the various subprocesses, 
we follow Ref. Jh}. 

Our kinematic cuts, imposed on the final state pho- 
tons, reflect the optimal detection capabilities of the AT- 
LAS detector j§: 

• p7p > 25 GeV, for the transverse momentum of each 
photon, 

• |y 7 | < 2.5, for the rapidity of each photon, and 

• p\j(j>\ +Pt) < 0.7, to suppress the fragmentation 
contribution, where p\ is the transverse momentum 
of the photon with the higher px value. 

Additionally, we restrict the invariant mass of the photon 
pair in the relevant region for the light Higgs production: 
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FIG. 2. Invariant mass distributions of photon pairs at 
the LHC, calculated using ResBos. The total curve (upper 
solid) is the sum of the 0(as) qq + qg — > 77V (dashed) and 
the 0(a%) gg + qg — > 77V (dotted) contributions. The lead- 
ing order curves for the contributions from O(ag) qq — » 77 
(dash-dotted) and 0{a 2 s ) gg — » 77 (lower solid) are also 
shown for comparison. 



80 GeV < Q < 150 GeV. We also apply a AR = 0.4 sepa- 
ration cut on the photons, but our results are insensitive 
to this cut (cf. Ref. @). 

In Fig. ^| we display the invariant mass distribution 
of photon pairs at the LHC in the inclusive process 
pp — > 77V, calculated using ResBos with the above cuts. 
We present separately the resummed contribution from 
the subprocesses 0(as) qq + qg — ^ 77V (dashed) which 
includes the fragmentation, and 0(a 3 s ) gg + qg — > 77V 
(dotted), as well as the total distribution (upper solid), 
which is the sum of these two curves. We also display, 
for comparison, the leading order contributions from the 
subprocesses O{a s ) qq — > 77 (dash-dotted) and 0(a 2 s ) 
gg — > 77 (lower solid). 

The normalization of the resummed cross-section is de- 
termined order-by-order by the coefficients C^ n \ with 
the coefficient derived from the NLO corrections. 
In Fig. H we used the exact for the subprocess 

qq+qg — > jjX, for which the complete NLO cross-section 
is known. For the subprocess gg + qg — > 77 A", in which 
the NLO virtual corrections have not yet been calculated, 
we used the approximate as described in the previ- 
ous Section. We found that the C(a|) cross-section is 
sensitive to . Namely, the resummed cross-section 
with the approximate included is about 1.9 times 
larger than the one with = 0. 

In the Q > 85 GeV region the contribution of the qq + 
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FIG. 3. Transverse momentum distributions of photon 
pairs at the LHC. The total resummed curve (upper solid) 
is the sum of the resummed qq + qg — » 77X (dashed), 
and the resummed gg + qg — > 77X (dotted) contributions. 
The fixed-order curves for the contributions from 0(as) 
99 + 19 ~ > 77^" (dash-dotted) and 0{a%) gg + qg — > 77X 
(middle solid) are also shown for comparison. The approxi- 
mation for the latter used in jlo| is also shown (lower solid). 

qg — > 77 subprocess is higher than that of the gg + 
qg — > 77V. From Fig. ^ we can also read the if -factors, 
which are defined as the ratios of the resummed to the 
leading order results. The A-factor for the qq + qg — > 
77V process (the ratio of the dashed and dash-dotted 
curves) is between 1.40 and 1.70 in the invariant mass 
range of interest. Similarly, the gg + qg — ► 77 V K- 
factor (the ratio of the dotted and lower solid curves) 
is between 1.45 and 1.75. This results in an overall K- 
factor of 1.4 to 1.7. These A-factors are about the same 
as the NLO/LO A"- factors in the fixed order perturbative 
calculations 

In Fig. we plot the transverse momentum distribu- 
tion of photon pairs at the LHC. In addition to the total 
resummed result, we give the resummed and fixed-order 
calculations separately for both the 0(as) qq + qg —> 
77A and the 0(a|) gg + qg — ► 77A subprocesses. In 
both cases the resummed results deviate substantially 
from the fixed order predictions in the < Qt < 100 
GeV region. At Qt — 30 GeV the resummed curves are 
higher by about 30 and 50 percent for the qq + qg — * 77 A 
and gg + qg — > 77A subprocesses, respectively. As a 
result the total resummed curve exceeds the total fixed- 
order prediction by almost 40 percent at Qt = 30 GeV. 
This is the Qt region where the kinematic cuts are ap- 
plied in order to optimize the statistical significance of 
the signal in the Higgs plus jet mode. Thus, the use of 
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FIG. 4. Invariant mass distributions of photon pairs 
at the LHC, with the cut Q T > 30 GeV. The total re- 
summed curve (upper solid) is the sum of the resummed 
99 + 19 — * 77^ (dashed) , and the resummed gg + qg —* 77 A 
(dotted) contributions. The fixed-order curves for the contri- 
butions from 0(ots) 11 + qg — » 77X (dash-dotted) and 0(a%) 
99 + 19 ~~* 77^f (lower solid) are also shown for comparison. 

the resummed prediction is necessary to extract a reli- 
able statistical significance, and also to make a correct 
determination of the Higgs production cross section, in 
the presence of kinematic cuts. 

Fig. H also shows that if the photon pair is constrained 
to be in the mid- to high-Qr region the contribution of 
the gg + qg — > 77V subprocess is small. At Qt = 40 
GeV, for example, the gg initial state accounts for less 
than 30 percent of the total cross section. In the ResBos 
program, the gg + qg — > 77V rate is predicted purely by 
the resummed calculation and does not cross over into the 
fixed-order 0(a%) calculation until after about Qt = 100 
GeV, at which point this rate is negligible. For reference 
we also show the approximate fixed order curve used in 
JToj l . The curve calculated using the exact matrix element 
significantly exceeds the approximation in the high Qt 
region. 

To illustrate the higher order effects on the invariant 
mass distribution in the presence of a Qt cut, in Fig. ^ 
we plot Q of the photon pair while restricting Qt > 30 
GeV. Due to the different shape of the Qt distributions 
the cut offsets the fixed-order and resummed rates, as 
explained in Ref. |ll[]. The effect is larger for the gg 
channel, since there the resummed and fixed-order Qt 
distributions deviate more, signaling higher corrections 
to the gg process than to the qq process. In the presence 
of kinematic cuts the difference in the Q distribution, 
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between the fixed order and the resummed calculations, 
can be as high as 50 percent. The Qt > 30 GeV cut also 
suppresses the gg channel, decreasing the uncertainty of 
the total prediction. 

IV. CONCLUSIONS 

In this paper we present the resummed calculation of 
the pp — > yyX distributions including the exact fixed or- 
der gg — > 773 contribution. Combining the known fixed 
order QCD corrections and the most important logarith- 
mic terms of the higher order corrections, the resummed 
cross section provides a reliable prediction for the inclu- 
sive diphoton invariant mass and transverse momentum 
(Qt) distributions. With a Qt cut the least reliable 
gg — ► 77A component can be suppressed, and the pre- 
diction further improved. 

Given our results, we propose the search for the Higgs 
boson in the inclusive diphoton mode with a cut on the 
transverse momentum of the photon pair. This measure- 
ment can be done without the requirement of a tagged 
jet, which is necessary in the 77 jet mode. Therefore, it 
is independent of the jet algorithm used, it can be per- 
formed more precisely experimentally, and it can be pre- 
dicted more reliably from a resummed calculation such 
as presented here. 

While finishing this paper we became aware of a similar 
work, in which the authors extract the 3g2"/ amplitude 
from the 5g amplitude fl6| . Our fixed order analytical 
and numerical results agree with the results of that paper. 
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